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Abstract 

The curvature coordinates T, R of a Schwarzschild spacetime are turned 
into canonical coordinates T(r),R(r) on the phase space of spherically 
symmetric black holes. The entire dynamical content of the Hamilto- 
nian theory is reduced to the constraints requiring that the momenta 
Pr(r), Pr(t) vanish. What remains is a conjugate pair of canonical vari- 
ables m and p whose values are the same on every embedding. The 
coordinate m is the Schwarzschild mass, and the momentum p the dif- 
ference of parametrization times at right and left infinities. The Dirac 
constraint quantization in the new representation leads to the state func- 
tional \&(ra;T, R] = ^(m) which describes an unchanging superposition 
of black holes with different masses. The new canonical variables may be 
employed in the study of collapsing matter systems. 



1 Introduction 

1.1 Context 

General relativity was cast into canonical form by Dirac and by Arnowitt, 
Deser and Misner (ADM) Inside asymptotic regions, Hamiltonian dynamics 
is entirely generated by constraints. Their imposition as operator restrictions on 
the states yields the Wheeler-DeWitt equation H , 0] • DcWitt realized that by 
freezing all but few degrees of freedom of a cosmological model by symmetry, one 
can obtain exactly soluble models of quantum gravity |Q . Misner and his school 
turned this idea of minisuperspace quantization H into a systematic exploration 
of quantum cosmology Q , |7| . 
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Minisuperspace techniques were extended to midisuperspace quantization of 
infinitely-dimensional models. The first system treated in this manner was the 
cylindrical gravitational wave ||. ft became clear that the Wheeler-DeWitt 
equation is often unwieldy and difficult to interpret. For cylindrical waves, the 
switch to an extrinsic time representation changed the Wheeler-DeWitt equation 
into a functional time Schrodinger equation. Gravity assumed the form of a 
parametrized field theory |Q, ||. 

Things did not work that way for other infinitely-dimensional systems. The 
most important of these is the gravitational collapse of a spherically symmetric 
distribution of matter. Berger, Chitre, Moncrief and Nutku (BCMN) set this 
problem in the Dirac-ADM midisuperspace formalism. In their classic paper 
| jio| they studied a spherically symmetric massless scalar field coupled to grav- 
ity. They did not succeed in finding an extrinsic time representation. Instead, 
they reduced the action to a privileged foliation characterized by the vanishing 
'radial' momentum. Their reduced Hamiltonian did not quite reproduce the 
field equations. This was found and corrected by Unruh |llj . 

The BCMN model opens a canonical route to the study of Hawking's radi- 



ation 12 , [|13|. A standard semiclassical analysis of the Hawking effect starts 
with a black hole being formed by the gravitational collapse of classical matter. 
One studies a field that propagates on this background. The modes that disap- 
pear below the horizon are averaged out, and the thermal radiation escaping to 
infinity is described by a density operator. The Hawking radiation leads to the 
evaporation of the black hole. No general agreement has been reached on what 
is the final state of this process. The black hole may evaporate completely, or 
leave a remnant. If it evaporates completely, the question remains what happens 
to the information which got initially trapped below its horizon. 

Midisuperspace canonical approach has two potential advantages over the 
standard description. First, it goes beyond the semiclassical approximation. 
Second, unless one encounters a Cauchy horizon, all information is registered 
in the canonical data on a Cauchy hypersurface. One can study what happens 
inside black holes and how they approach the singularity. However, to make 
use of this advantage, one must insist that the foliation covers all available 
spacetime, and that time evolution is not artificially arrested. The BCMN 
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slicing does not meet this condition. A careful study of the BCMN model in 
its relation to Hawking's radiation was undertaken by Hajfcek - [[l6|. He 
generalized the model to other spherically symmetric fields, and paid special 
attention to the properties of apparent horizon and choice of slicing. 

One can easily see that the BCMN slicing is problematic already for primor- 
dial Schwarzschild black holes. In vacuo, the BCMN slices coincide with those 
of constant Killing time T. They cover only the static regions of the Kruskal 
diagram, and never penetrate the horizon. Canonical treatment of a complete 
Schwarzschild spacetime was attempted by Lund jl7| . To get below the horizon, 
Lund used the Lamaitre slices, or the slices of constant R < 2M. He did not 
succeed in covering the whole Kruskal diagram by a single foliation, or relate 
the state of the Schwarzschild black hole on T = const slices to its state on the 
Lamaitre slices or the 2M > R = const slices. The best solution would have 
been to work in the functional time representation, but Lund presented a proof 
that the extrinsic time representation does not exist for vacuum Schwarzschild 
black holes. 

Interest in primordial black holes has been revived by a surge of activ- 
ity on quantization of dilatonic black holes (see recent reviews by Gidddings 
[ fl8| , and by Harvey and Strominger [|l9| ) . Starting from this program, Gegen- 
berg, Kunstatter and Louis-Martinez |^0| , plj] discussed canonical quantization 
of Schwarzschild black holes within conformally invariant formulation of Ein- 
stein's theory |22]. Along different lines, Thiemann and Kastrup [^3| - [^6| 
discussed canonical quantization of Schwarzschild and Riesner-Nordstrom black 
holes mostly, though not entirely, within Ashtekar's canonical formalism. 

Keeping this background in mind, let us state the goals and results of this 
paper. 

1.2 Results 

We cast the classical and quantum dynamics of primordial black holes into geo- 
metrically transparent and explicitly soluble form by choosing a natural canon- 
ical chart on the phase space. Our method amounts to finding a functional 
extrinsic time representation analogous to one which exists for cylindrical gravi- 
tational waves. Lund's no-go theorem is transcended because the representation 
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does not satisfy its unnecessarily strong premises. 

Geometrodynamics of Schwarzschild black holes is governed by the Dirac- 
ADM action restricted to spacetimes with spherical symmetry. The spatial 
metric g a b(z c ) on a symmetric hypersurface is entirely characterized by two 
functions, A(r) and R(r), of a radial coordinate r. In canonical formalism, these 
are accompanied by the conjugate momenta P\(r) and Pr{t). We choose phase 
space variables which have an immediate geometric meaning. The hypersurface 
action yields the familiar Hamiltonian and momentum constraints equivalent to 



those derived by BCMN (10 



Proper understanding of the canonical formalism requires a careful handling 
of boundary conditions. One must specify how fast the canonical data and the 
lapse-shift multipliers fall at infinities. The requirement that A can be freely 
varied within its falloff class mandates the addition of the ADM boundary term 
to the hypersurface action. Unfortunately, the introduction of this term prevents 
one from freely varying the lapse function at infinities. This defect is amended 
by parametrization, which makes the action dependent on two more variables, 
namely, on proper times measured by static clocks at infinities. All the variables 
in the parametrized action can be freely varied, and their interconnection at 
infinities acquires the status of natural boundary conditions. This version of 
the canonical action principle is vital for our treatment of black hole dynamics. 

The crux of our approach is the introduction of the Killing time T(r) as 
a canonical coordinate. The way in which T(r) enters the formalism is rather 
delicate. To begin with, one restricts attention to the space of solutions. In 
a given Schwarzschild spacetime, one can specify a hypersurface by giving the 
familiar curvature coordinates of its points as functions of a radial label r. 
The Schwarzschild geometry induces on the hypersurface the canonical data 
which satisfy the constraints. From those data, one can locally determine the 
Schwarzschild mass M of the embedding spacetime, and the rate T'(r) at which 
the Killing time changes along the hypersurface. Though T(r) becomes infinite 
on the horizon, its change across the horizon can be consistently inferred from 
smooth canonical data. 

Next, one forgets how the expressions for M{r) and —T'(r) were obtained, 
and turns them into definitions of two new sets of dynamical variables on the 
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phase space. At this stage, the canonical data no longer need to satisfy the 
constraints, and the mass function M(r) can in principle depend on r. The 
remarkable feature of the new variables is that they form a canonical pair: 
Pm{t) — —T'(r). Moreover, by retaining the curvature coordinate, R(r) <— > 
R(r) = R{r), but modifying its conjugate momentum, Pr(t) i— ► Pr(t*), one can 
complete M(r) and —T'(r) into a new canonical chart 

M(r), P M (r) = -T'(r) ; R(r) = R(r), P R (r) . (1) 

We exhibit the generating functional of the canonical transformation from the 
old to the new canonical variables. Because T'(r) has an infinite jump on the 
horizon, the canonical transformation has there a singularity. 

We can now reimpose the constraints, and evolve the variables by Hamilton 
equations. Under these circumstances, the mass function becomes a position- 
independent constant of motion. Inversely, one can show that the Hamiltonian 
and momentum constraints can be replaced by a much simpler set of conditions 

M'(r) = 0, P R (r)=0 (2) 

on the new canonical variables. 

The canonical structure we have derived holds only when the lapse func- 
tion in the ADM boundary action is considered as a fixed function of the 
label time. After parametrization, the action becomes dependent on a pair 
of proper times r± at infinities, and it no longer has a canonical form. It is 
surprising that without adding any more variables to the parametrized space 
r+, T-, M(r), Pm(t), R(0, Pk( r )i one can introduce on it a canonical chart 

m,p; T(r),P T (r); R(r),P R (r). (3) 

The final canonical variables (j|) have a simple physical meaning. The canonical 
pair R(r),P R (r) remains unchanged. The new canonical coordinate T(r) is the 
Killing time, and its conjugate momentum Pt{t) is the mass density —M'(r). 
The curvature coordinates T, R in spacetime are thereby turned into canonical 
coordinates T(r), R(r) on the phase space. Simultaneously, the constraints are 
transformed into the statement that the momenta Pt(t) and Pr(t) canonically 
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conjugate to the embedding variables T{r) and R(r) vanish: 

P T (r) = 0, P R (r)=0. (4) 

The Hamiltonian is a linear combination of these constraints. The true ADM 
Hamiltonian disappeared in the transformation process to the final canonical 
chart (0) . This phenomenon can be understood as a result of time-dependent 
canonical transformation. 

Because the Hamiltonian weakly vanishes, the remaining canonical vari- 
ables m and p are constants of motion. The canonical coordinate m is the 
Schwarzschild mass of the black hole. The meaning of p is more esoteric: p char- 
acterizes the difference between the parametrization times at the left and right 
infinities. The comparison of these two times is made possible by connecting 
the infinities by hypersurfaces of constant Killing time. Once set, this difference 
is preserved because the two parametrization clocks run at the same rate. This 
explains why p is a constant of motion. While it has always been suspected that 
the variable conjugate to the Schwarzschild mass is some sort of proper time, 
the correct interpretation of this quantity eluded previous investigators, likely 
because they did not pay enough attention to the role of parametrization, and 
to how the two asymptotic regions are connected with each other. 

The parametrized action can be reduced to true dynamical degrees of free- 
dom, and the dynamics generated by the unparametrized action can be com- 
pared to that we have just described. It transpires that our parametrized view- 
point corresponds to performing a time-dependent canonical transformation to 
'initial data'. 

With the new polarization (U) of the phase space, the Dirac constraint quan- 
tization of primordial black holes becomes straightforward. The state func- 
tional ^(m; T, R] of the system is reduced by the constraints to an embedding- 
independent function \l/(m) of the mass parameter m. Such a state function 
describes a superposition of primordial black holes of different masses. Once 
prepared, it stays the same on every hypersurface T(r), R(r). 

The curvature coordinates are ill-behaved on the horizon. However, once 
they are constructed from the canonical data, one can easily transform them 
into Kruskal coordinates by a further canonical transformation. 
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Primordial black holes, despite all the care needed for their proper canoni- 
cal treatment, are dynamically trivial. The true interest of the new canonical 
variables lies in the possibilities which they open in the study of gravitational 
collapse of matter. These questions are now being pursued in collaboration with 
Dr. Petr Hajfcek and Dr. Joseph Romano. 

2 Schwarzschild Solution 

2.1 Spacetime Description 

The recurrent theme of this paper is that canonical formalism should be guided 
by spacetime intuition. We thus start by summarizing what is known about 
spherically symmetric solutions of vacuum Einstein equations. 

Any such solution is locally isometric to the Schwarzschild line clement 

ds 2 =-F(R) dT 2 + F~ 1 (R) dR 2 + R 2 dtt 2 (5) 

written in the curvature coordinates (T, R) . Here, 

dn 2 = dO 2 + sin 2 6 d(f> 2 (6) 

is the line element on the unit sphere. We are using natural units in which the 
Newton constant of gravitation G and the speed of light c are put equal to one: 
G = 1 = c. The coefficient F(R) has the form 

F(R) = l-2M/R, (7) 

where M is a constant. The curvature coordinate R is invariantly defined by the 
requirement that AttR 2 be the area of the 2-spheres S 2 : T = const, R = const 
which are the transitivity surfaces of the rotation group. The vector field d/dT 
is a Killing vector field of the metric (jE)). It is orthogonal to the hypersurfaces 
T = const of the Killing time T. 

For M = the spacetime is flat. Solutions with M > describe black holes, 
solutions with M < correspond to naked singularities. We limit ourselves to 
solutions with M > 0. 

As R — > oo, the line element (||) becomes asymptotically flat. At R = 
2M, the solution (0) runs into a coordinate singularity. The maximal analytic 
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extension of (p|) from the region R > 2M > across R — 2M describes a 
primordial black hole. The complete spacetime M. is represented by the familiar 
Kruskal diagram |27| , [^8| . It is covered by four patches of curvature coordinates 
which meet at the horizon F(R) = 0. Two regions have R > 2M, and the Killing 
field d/dT in them is timelike. We call them the right and left static regions, I 
and III. Two other regions have R < 2M, and the Killing field d/dT in them is 
spacelike. We call them the past and future dynamical regions, IV and II. The 
past dynamical region begins and the future dynamical region ends in a true 
curvature singularity at R = 0. 

The total Schwarzschild spacetime can be covered by a single patch of 
Kruskal coordinates U and V. The lines U — const are radial rightgoing null 
geodesies /*, the lines V — const are radial leftgoing null geodesies \. Both U 
and V grow from past to future. The horizon is transformed to the lines U = 
and V = 0. Their intersection is the bifurcation point. The four regions covered 
by curvature coordinates are: 

right static region I: R > 2M : U < 0, V > , 

left static region III: R > 2M : U > 0, V < , 

future dynamical region II: R < 2M : U > 0, V > , 

past dynamical region IV: R < 2M : U < 0, V < . 

The Kruskal coordinates U and V are mapped into curvature coordinates 
T and R by a two-to-one transformation. Anticipating the steps we shall need 
later in the canonical formalism, we build this transformation in several steps. 
The Kruskal coordinates are dimensionless. We thus first scale the curvature 
coordinates T and R into dimensionless coordinates T and R: 

T=—,R=—. (8) 

2M 2M v ' 

Curvature coordinates are related to Kruskal coordinates by 

UV = W{R) := (1 - R) exp (R) , - = sign(l - R) exp (f) . (9) 

The T = const hypersurfaces appear in the Kruskal diagram as straight lines 
passing through the bifurcation point. The rightgoing branch /" of the horizon 
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is labeled by T = oo, the lcftgoing branch \ by T = — oo. The R — const lines 
are hyperbolas asymptotic to the horizon. 

As R increases from R = (singularity) to R = oo, the function W(R) 
monotonically decreases from 1 to — oo. Therefore, it has an inverse, 72.(17 V). 
This enables us to solve (^|) for the curvature coordinates: 

f = ln|V| -ln|Z7|, i? = ft(UV) . (10) 

We see that the inversion U t— > — C/,V^i— > — V leaves T and 7? unchanged. Two 
points in the Kruskal plane are labeled by the same curvature coordinates. 



It is useful to describe the mapping (10) in a slightly weaker form. Introduce 
the tortoise coordinate 

R* = R + ln\ 1 -R\ (11) 
and combine T and R* into null coordinates 



U = T-R*, V = T + R*. (12) 

By multiplying and dividing the two equations (0) we obtain 

U 2 =exp(-U), V 2 = cxp (V) . (13) 

The Schwarzschild line element expressed in the Kruskal coordinates is every- 
where regular, except at the initial and final curvature singularities at UV = 1. 

2.2 Geometrodynamical Description 

Geometrodynamics views a given spacetime as a dynamical evolution of a three- 
geometry. Let us apply this viewpoint to the Schwarzschild solution (0). Take 
an arbitrary spherically symmetric spacelike hypersurface and let the spacetime 
metric induce on it the spatial geometry g. Different hypersurfaces carry differ- 
ent induced metrics. As we evolve the spacetime metric along a foliation which 
covers the entire Kruskal diagram, we build the Schwarzschild solution. 

The simplest evolution is obtained along a one-parameter family of hyper- 
surfaces T = const cutting across the regions I and III of the Kruskal diagram. 
The geometry induced on all of these hypersurfaces is the same: it is a wormhole 
geometry known under the name of the Einstein-Rosen bridge. As we change T 
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from — oo to oo, the geometry does not change at all. It would be more appropri- 
ate to speak in this case about geometrostatics rather than geometrodynamics. 
This reflects the fact that d/dT is a Killing vector field, and the T = const 
hypersurfaces are Lie-propagated by it. This evolution has a serious defect p9| : 
The spacelike hypersurfaces T — const do not cover the entire Kruskal diagram, 
but only the static regions I and III. The progress of time is arrested at the 
bifurcation point through which all the hypersurfaces pass. The hypersurfaces 
thus do not form a foliation. Moreover, as the hypersurfaces proceed from past 
to future in region I, they recede from future to past in region III. The evolution 
does not proceed everywhere from past to future. 

The dynamical regions II and IV which were not covered by the spacelike 
hypersurfaces T = const can be covered by another simple family of spacelike 
hypersurfaces, namely, R = const < 2M. Their geometry is again a wormhole 
geometry, but this time it describes a homogeneous cylinder S 2 xJR. As R pro- 
gresses from to 2M in region IV, the cylinder opens up from the line singularity, 
its circumference grows larger and larger while its length per unit T shrinks, 
and finally degenerates into a disk of circumference AttM as R approaches the 
horizon R = 2M. In region II, as R decreases from 2M to 0, the whole process 
is reversed. Spatial geometry is dynamical in regions IV and II. The described 
evolution is locally isomorphic to the dynamics of the Kantowski-Sachs universe 
[ p0| . Unfortunately, the foliations < R = const < 2M are asymptotically 
null at infinities, they approach the horizon for R — > 2M , and again, they do 
not cover the whole Kruskal diagram. What they miss are exactly the static 
regions. 

These shortcomings set our task. We want to study the spatial geometry 
on a hypersurface which is spacelike and cuts the Kruskal diagram all the way 
through: it starts at left infinity, goes through the static region III, crosses the 
horizon into a dynamical region, traverses it until it again reaches the horizon, 
crosses the horizon to the static region I, and continues to right infinity. At 
infinities, such a hypersurface should be asymptotically spacelike, approaching 
some static hypersurface T = T_ = const at left infinity, and in general some 
other static hypersurface, T = T + = const, at right infinity. One can cover the 
whole Kruskal diagram by a foliation of such hypersurfaces. Even better, one 
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can admit all of them at once, and work in many-fingered time formalism. 

3 Canonical Formalism for Spherically Symmet- 
ric Spacetimes 

The geometrodynamical approach does not start from the known Schwarzschild 
solution: it generates the Schwarzschild solution by evolving a spherically sym- 
metric geometry. The evolution is governed by the Dirac-ADM action. In this 
section, we introduce the Dirac-ADM action, and carefully discuss the necessary 
boundary conditions. 

3.1 Hypersurface Lagrangian 

Take a spherically symmetric three-dimensional Riemannian space (E, g) and 
adapt the coordinates x a of its points x € E to the symmetry: x a = (r,8,<j)). 
The line element da on E is completely characterized by two functions, A(r) 
and R(r), of the radial label r: 

da 2 = A 2 (r) dr 2 + R 2 (r) dil 2 . (14) 

Again, dfl is the line element on the unit sphere. 

Primordial black holes have the topology E = BxS 2 , and r G 1R thus ranges 
from — oo to oo. The coefficients A(r) and R(r) cannot vanish because the line 
element must be regular. We take both of them to be positive, A(r) > and 
R{r) > . Then, R(r) is the curvature radius of the 2-sphere r = const , and 
da = A(r)dr is the radial line element oriented from the left infinity to right 
infinity. 

Under transformations of r, R(r) behaves as a scalar, and A(r) as a scalar 
density. This will simplify the momentum constraint. Keeping this goal in 
mind, we have avoided the usual exponential form of the metric coefficients. It 
is important to keep track of the density character of the fundamental canonical 
variables. We shall always denote those canonical coordinates which are scalars 
by Latin letters, and those which are scalar densities by Greek letters. 

The line element (|l4]) leads to the curvature scalar 

R[g] = -AA^R^R" + 4A" 3 i?" 1 A'i?' - 2A" 2 ir 2 i?' 2 + 2R~ 2 . (15) 



11 



Let us foliate a spherically symmetric spacetime M by spherically symmetric 
leaves £, and label the leaves by a time parameter t € IR. The metric coefficients 
A and R then depend not only on r, but also on t. The leaves of the foliation 
are related by the familiar lapse function N and the shift vector N a . Due to the 
spherical symmetry, only the radial component N r of the shift vector survives, 
and both N(t, r) and N r (t, r) depend solely on the t, r variables. 

The extrinsic curvature K a b of the leaves is given by the rate of change g a b 
of the metric with the label time t : 



K ab = 9ab + #(o|5)) • 



For the spherically symmetric line element (114), K a b is diagonal, with 



Kr 



-N~ 1 a(A- (AiV r )') 



N~ R[ R — R N 1 



sin 2 6 Kg 9 . 



The vacuum dynamics of the metric field follows from the ADM action 



(16) 

(17) 
(18) 

i 

Sz[g,N,N a } = J dt J d 3 xL^ (19) 

whose Lagrangian is 

is = (lG^r 1 N\g\^ 2 (K ab K ab -K 2 + R[g}) . (20) 

(In natural units G = 1 = c, the Einstein constant k = 87rG/c 4 reduces to 8ir.) 
In a spherically symmetric spacetime (14), jl7|), ( Eg ) 

N\g\ 1/2 {K ab K ab -K 2 ) = -27V- 1 sin 61 

x (2(- A+(AN r )')(- R + R'N r )R+ (- R + R'N r fA) . (21) 

Integration over 8 and <\> gives the ADM action of a Schwarzschild black hole: 

//■oo 
dt / dr 
J —OO 



N~ 



■ N 



R(-A+ (AN r )') (-R + R'N r ) + ^A( -R+ R' N r f 



-A^RR" + A- 2 RR'A' - ^A^R' 2 + ^A 



(22) 



We shall discuss the appropriate boundary terms after passing to the Hamil- 
tonian formalism. 



12 



3.2 Canonical Form of the Action 

By differentiating the ADM action with respect to the velocities A and R we 
obtain the momenta 

Pa = -N- 1 R(h- R'N r ^j , (23) 

P R = -N- 1 (A(R - R'N r ) + R(A - (AAH')) • (24) 

Throughout this paper, we denote those canonical coordinates which are spatial 
scalars by Latin letters, and those which are spatial densities by Greek letters. 
The conjugate momenta always carry complementary weights. Therefore, the 
momentum P R conjugate to the scalar R is a density, while the momentum Pa 
conjugate to the density A is a scalar. 

Equations (^3|) and ( p4| ) can be inverted for the velocities, 

A = -NR- 2 (RP R - AP A ) + (AN r )' , (25) 
it = -N R- 1 P A + R'N r . (26) 

They allow us to write the extrinsic curvature as a function of the canonical 
momenta: 

K rr = AR- 2 (RP R -AP A ), K ee = P A . (27) 

By symmetry, the curvature of a normal section of £ attains its extremal 
values Ki, K 2 , K 3 (called principal curvatures) for those sections which are ei- 
ther tangential || or normal _L to the 2-spheres r = const. The last equation 
enables us to express these principal curvatures K\ = K\\ and K 2 = K§ = K± 
in terms of the momenta: 

K\\ = K = A- 1 R-' 2 {RP R - AP A ) , 

Kx_ = K e g =Kl = R- 2 P A . (28) 

Inversely, 

P A = R 2 K ± , P R = RA {K\\ + Kj_) . (29) 

This endows the canonical momenta with an invariant geometric meaning. Note 
that Pa is proportional to but P R is proportional to the sum of Ku and 
K\_ , not to if || itself. This sum also differs from the mean curvature 

K :=K 1 +K 2 + K 3 =K2 = K }] + 2K ± . (30) 
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The action (|22j) can be cast into the canonical form by the Legendre dual 
transformation: 

Ss [A, R, Pa, Pr; N, N r ] = f dt f dr (P A A + P R R - NH - N r H r ) . (31) 
By this process, we obtain the super-Hamiltonian 
H:= -R- x P R P K +^R- 2 KPl 

+ A~ l RR" - A- 2 RR'A' + ^A^R' 2 - ^A (32) 

and supermomentum 

H r :=P R R'- AP A . (33) 

The form of H r is dictated by the requirement that it generate DifflR of the 
scalars R and Pa, and of the scalar densities Pr and A. The minus sign in ( |33"| ) 
is due to the fact that the momentum Pa is a scalar, and the coordinate A a 
scalar density, rather than the other way around. 

The expressions ([32]) and ( |33| ) can be obtained from those derived by BCMN 
[ flO| by a point transformation. 

3.3 Falloff of the Canonical Variables 

So far, we paid no attention to the behavior of the canonical variables A, R and 
Pa, Pr at infinity. The importance of the falloff conditions was pointed out by 
Regge and Teitelboim fl3l[ , |32]| , and their form refined by many authors. We 
shall follow the treatment of Beig and O'Murchadha |33[] . 

Primordial black holes have two spatial infinities rather than just one. We 
shall formulate the falloff conditions at right infinity, and then state what the 
corresponding conditions are at left infinity. 

Let x a be a global system of coordinates on S which is asymptotically Carte- 
sian. Such a system is related to the spherical system of coordinates r, 9, <f> by 
the standard flat space formulae. At r — > oo, the metric g a b and the conjugate 
momentum p ab are required to fall off as 

9ab(x c ) = <U + r- 1 /i ab (n c ) + 0°°(r-( 1+e >), (34) 
p ab (x c ) = r- 2 k ab (n c ) + oo (r- (2+e '>) . (35) 
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Here, n a :—x a /r, and f(x a ) = 0°°(r n ) means that / falls off like r ™, / )Q 
like r - ( - n+1 \ and so on for all higher spatial derivatives. The leading term in 
9ab — Sab is of the order r -1 , and the leading term in p ab is one order higher, 
r~ 2 . The coefficients h ab and k ab are required to be smooth functions on S 2 , 
and to have opposite parities: h ab (n c ) is to be even, and k ab (n c ) is to be odd, 
i.e. 

h ab {-n c ) = h ab (n c ) , k ab (-n c ) = -k ab (n c ) . (36) 

Together with the canonical data, the lapse function and the shift vector are 
assumed to behave as 

N(x c ) = N + (n c ) + 0°°(r- e ), (37) 
N a {x c ) = Nl{n c ) + 0°°{r- c ) (38) 

at infinity. 

Because dr = n a dx a and dVL 2 = r~ 2 (5 ab — n a nb) dx a dx b , the spherically 
symmetric metric ( p"4| ) can easily be transformed into Cartesian coordinates: 

g ab = A 2 n a n b + f (S ab - n a n b ) , (39) 

ib a —2^a^b ( \ I cab 



g ab = A- 2 n a n b + (-J (S ab - n a n b ) . (40) 

The general falloff conditions (Q) then determine the behavior of the metric 
coefficients R and A at infinity: 

A(t,r) = 1 + M+i^r- 1 + 0°°(r-^) , (41) 

and 

R{t,r) = r + p+{t)+O x (r- e ). (42) 

Here, due to spherical symmetry, M + (t) and p+ (t) cannot depend on the angles 
n c , but they can still depend on t. Of course, M + is the Schwarzschild mass 
as observed at right infinity. The general falloff conditions allow R to differ by 
the amount p+{t) from r at infinity . Because we want r to coincide with R at 
infinity, we impose stronger falloff conditions on R by putting p + = : 

R{r)=r + O x {r- e ). (43) 
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Consistency then demands that the shift vector must also asymptotically vanish, 



iV; = 0: 



N r {r) =0°°(r- £ ). (44) 



Unlike N r , the lapse function cannot vanish at infinity. If it did, the time there 
would stand still. Instead, N(t, r) assumes at infinity an angle-independent 
value N+(t) : 

N(r)=N + (t)+0°°(r- e ). (45) 

Our next task is to determine the falloff of the canonical momenta. The 
second fundamental form K ab dx a dx b can be transformed into Cartesian coordi- 
nates by the same procedure as the metric. By using (p7|) we get 



K rr dr 2 + K ge d9 2 + K^dcf 2 = 

AiT 2 (RP R - AP A ) dr 2 + P A dfl 2 = 

(AR- 2 (RP R - AP A )n a n b + r- 2 P A (5 ab - n a n b )^j dx a dx b . (46) 

From here we read the extrinsic curvature and, by referring back to the metric 
(^0|), find the canonical momentum 

p ab :=\g\ 1/2 {Kg ab ~ K ab ) 

= 2r- 2 A- 1 P A n a n b + R- 1 P R (d ab - n a n b ) . (47) 

This expression for p ab is even in the angular variables. This means that k ab is 
even. The requirement that k ab be odd implies that k ab must vanish, i.e. 

P A =0°°(r- e ), P H = oo (r-( 1 + £ )). (48) 

The same considerations apply at r — > — oo . Taken together, they imply the 
falloff conditions 

A(t,r) = l + M ± (t)|r|- 1 +0 oo (|r|-( 1+e )), (49) 
R(t,r) = |r|+0°°(|rn, (50) 
P A (t,r) = 0°°(|r|- £ ), (51) 
P R (t,r) - 0°°(|r|-( 1+e )) (52) 
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for the canonical variables, and the falloff conditions 



N{t,r) = N±(t)+0°°{\r\- e ), (53) 
N r (t,r) = 0°°(\r\- e ) (54) 

for the Lagrange multipliers. 

Let us note again that R(r) approaches \r\ at the rate 0°°(|r|~ c ) as r — > ±00. 
After we reconstruct the Killing time T(r) from the canonical data, we shall 
prove a similar result for the approach of the t = const foliation to the T = const 
foliation: 

T(t,r) =T ± (t)+0°°(\r\-). (55) 
The falloffs of the canonical data ensure that the Liouville form 

/oo 
dr(P A A + P R R) (56) 
-00 

is well defined. They also imply that the super-Hamiltonian and supermomen- 
tum fall off as 

H = 00 (\r\- {1+e) ) 1 H r = 0°°(\r\- {1+e) ). (57) 

Equations ([33]), ( pi| ) and (|57]) ensure that the Hamiltonian is well defined. The 
canonical action Ss thus has a good meaning. 

The falloff conditions for spherically symmetric vacuum spacetimes may eas- 
ily be strengthened. The only necessary condition on how fast the canonical vari- 
ables fall down is that the Schwarzschild solution itself can be made to fall down 
that fast by an appropriate choice of the hypersurface and its parametrization 
r at infinities. By choosing the hypersurface to coincide with a T — t hypersur- 
face outside a spherical tube R = Rq = const , and by labeling it there by the 
curvature coordinate, R = \r\ , we achieve that 

A(t,r) = l + M ± (t)\r\- 1 , R(t,r) = \r\, 

P A (t,r)=0, P«(t,r) = ( ™ ' 

for the Schwarzschild solution. We are thus free to require that our phase space 
variables satisfy (pa) everywhere outside a compact region. 
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3.4 Boundary Terms 

Arnowitt, Deser and Misner complemented the hypersurface action Ss by a 
boundary action Son at infinity: 

S = + S as ■ (59) 

A primordial black hole has two infinities, and hence there are two boundary 
contributions 



S as = - / dt (N+(t)E+(t) + N-(t)E-(t)) . (60) 
Each of them is the product of the lapse function with the ADM energy 

E± = (167T)- 1 f dS a 5 bc (g ab . c - g bc , a ) (61) 

which is given by an integral over a two-sphere Sj_ at infinity. In the asymptot- 
ically Cartesian coordinates 

dS a = r 2 n a sin 9 d6d<p, \r\. a =n a , n a . b = \r\- 1 (6 ab - n a n b ) , (62) 

and the asymptotic form ( |34| ) of the metric yields 

E±= lim i|r|(A 2 -l) =M±(t). (63) 

The ADM energy of a black hole is its Schwarzschild mass. 
The total action of a primordial black hole takes the form 

S[A,P A ,R,P R ; N, N r ] = 

dt dr (P A A + P R R - NH - N r H r ) 

J J —OO 

dt (N+M+ + iV_Af_) . (64) 



Variations of the canonical variables should preserve the prescribed falloffs. 
In particular, the leading term in the variation of A is given by the variation of 
the Schwarzschild mass: 

5A± = 5M±\r\- 1 +0°°(\r\- 2 ). (65) 
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As emphasized by Regge and Teitelboim [^IJ , without the ADM boundary action 
this particular variation of the action would lead to an inconsistency. On the 
other hand, when the ADM boundary action is included, another inconsistency 
would arise if one allowed the variation of the lapse function at infinities. It is 
therefore important to treat N±(t) as prescribed functions of t. 

To see why this is so, let us identify those terms in the Hamiltonian density 
NH + N r H r whose variations lead to boundary terms. The variation of Pr does 
not yield any boundary term because there are no derivatives of Pr anywhere 
in the action. The variables Pa and R also do not cause any trouble, because 
the falloff conditions ensure that the boundary terms brought in by the varia- 
tion of Pa and R safely vanish. The sole troublemaker is the derivative term 
NRR'K~ 2 A! in —NH . Its variation with respect to A yields the boundary term 

N+(t)5M+(t) + N-(t)5M-(t) . (66) 

If there were no ADM boundary action, the variation of the hypersurface action 
Ss with respect to A would lead to the conclusion that N±(t) = , i.e. it would 
freeze the evolution at infinities. The ADM boundary action is designed so that 
its variation with respect to M±(t) exactly cancels the boundary term (p6|), and 
the unwanted conclusion does not follow. 

The ADM boundary action leads, however, to an inconvenient caveat. With- 
out it, the lapse and the shift functions in the hypersurface action Ss can be 
freely varied at infinities because the super-Hamiltonian H and the supermo- 
mentum H r asymptotically vanish. When the boundary action is included, the 
variation of the total action S = Ss + Sgs yields 

S N S = - J dt (M+(t)SN + (t) + M_ (t)5N-) . (67) 

If we allowed N(t, r) to be varied at infinities, we would get an unwanted 'natural 
boundary condition' M+ = = M_. This would exclude black hole solutions 
and leave only a flat spacetime. Therefore, in the variational principle ( |64| ) we 
must demand that the values N± (t) of the lapse function at infinities be some 
prescribed functions of t which cannot be varied. This means that the lapse 
function N in our variational principle has fixed ends. 
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3.5 Parametrization at Infinities 

The necessity of fixing the lapse function at infinities can be removed by paramet- 
rization. The lapse function is the rate of change of the proper time r with 
respect to the label time t in the direction normal to the foliation. Because 
N±(t) = 0, we can write 

N±(t) = ±r±(t) , (68) 

where r±(t) is the proper time measured on standard clocks moving along the 
r = const worldlines at infinities. By convention, we let the proper time at the 
left infinity decrease from the past to the future, to match the behavior of the 
Killing time T in the Kruskal diagram. This introduces the minus sign in ( |68|) 
at — oo. 

We now replace the lapse function in the ADM boundary action by the 
derivatives of t±, (|68|), and treat T±(t) as additional variables: 

S[A,P A ,R,P R ; N, N r ; r+,r_] = 

f dt dr (P A A + P R R - NH - N r H r ) 

J J — OO 

- J dt (M+t+ - M_f_) . (69) 

This rearrangement of the action is called the parametrization at infinities. No- 
tice that N± still appears in the hypersurface part 5*s of the action. 

The variables r±(t) in the action ( |6^ ) can be freely varied. The result of 
their variation is a valid equation, namely, the mass conservation 

M+(t) = = M-(t) . (70) 

The lapse function at infinities N±(t) can also be freely varied, because it now 
occurs only in the hypersurface action, and the super-Hamiltonian asymptoti- 
cally vanishes. 

The only remaining question is what happens under the variation of A. As 
before, the variation of the hypersurface action Ss gives the boundary term 
(p6[). On the other hand, the variation of the parametrized ADM boundary 
action in ( |69| ) now yields 

- T+(t)5M+(t) + T-(t)SM-(t) . (71) 
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Before parametrization, the variation of A at infinities, i.e. the variation of 
M±(t), produced merely an identity. After parametrization, the situation is 
different: the variation of M±(t) relates r±(t) to N±(t) by (|6|). These equations 
thereby follow as natural boundary conditions from the parametrized action 
principle (p9|). 

As we proceed, we shall at first pay little attention to the ADM boundary 
action. We shall meet it again in Section 

4 Reconstructing the Mass and Time from the 
Canonical Data 

4.1 Reconstruction Program 

In canonical gravity, we know the intrinsic metric and extrinsic curvature of 
a hypersurface, but we do not have any a priori knowledge about how the 
hypersurface is located in spacetime. Suppose that we are given the canoni- 
cal data A, R, P\, Pr on a spherically symmetric hypersurface cutting across a 
Schwarzschild black hole. Can we tell from those data the mass of the black 
hole? And can we infer from them how the hypersurface is drawn in the Kruskal 
diagram? 

Disregarding some subtleties concerning the anchoring of the hypersurface at 
infinities, the answer to both of these questions is yes. To arrive at the answers, 
we start from the knowledge that the hypersurface must ultimately be embedded 
in a spacetime endowed with the line element (||). Let the hypersurface be a 
leaf of a foliation 

T = T(t,r), R = R(t,r). (72) 

The line element (||) induces an intrinsic metric and extrinsic curvature on 
the hypersurface. By comparing these quantities with the canonical data, we 
connect the spacetime formalism with the canonical formalism. This connection 
enables us to identify the Schwarzschild mass M and the embedding (|72| ) from 
the canonical data. 

We can expect discontinuities at the horizon where different patches of the 
curvature coordinates T and R meet each other. Indeed, T becomes infinite at 
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the horizon. We shall see, however, that the transition of T across the horizon 
is well under our control, and can be predicted from smooth canonical data. We 
can then easily pass from the curvature coordinates T, R to the Kruskal coordi- 
nates U, V which are continuous across the horizon. The direct reconstruction 
of the Kruskal coordinates would be much more cumbersome. 

To implement our program, we substitute the foliation (|7^) into the Schwarz- 
schild line element (|^) and get 

ds 2 = -(FT 2 - F^R 2 ) dt 2 

+ 2( - FT'f + F~ l R'R) dtdr 

+ (- FT' 2 + F^R' 2 ) dr 2 + R 2 dn 2 . (73) 

By comparing this result with the ADM form of the line element 

ds 2 = -{N 2 -K 2 {N r f)dt 2 
+ 2k 2 N r dtdr 

+ A 2 dr 2 + R 2 dn 2 , (74) 

we obtain a set of three equations, 

A 2 = -FT' 2 + F- 1 R' 2 , (75) 
A 2 N r = -FT'f + F- 1 R'R, (76) 
N 2 -A 2 (N r ) 2 = FT 2 -F- X R 2 . (77) 

The first two equations can be solved for N r , 

—FT'f + F~ 1 R'R 

_ FT ,2 +F -l R ,2 ■ 

This solution N r , together with A of equation J75|), can be substituted into the 
remaining equation (f77|), which can then be solved for the lapse: 

R'f-T'R 

N = ■ . 79 

V-FT' 2 + F~ l R' 2 K ' 

The lapse function TV should be positive. The transition from (^) to ([79]) 
requires taking a square root. We must check that the square root we have 
taken is positive. 
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First of all, the denominator of ( |79[ ) is real and positive: it is equal to A. We 
shall check that the numerator of J79| ) is positive separately in each region of the 
Kruskal diagram. Because the lapse function is a spatial scalar, we first choose 
in each region an appropriate radial label r : r = R in region I, r = T in region 
II, r = —R in region III, and r = —T in region IV. Under these choices, the 
numerator of ( f79[) becomes T in region I, —R in region II, — T in region III, and 
R in region IV. With the label time going to the future, all these expressions 
are positive. 

We now substitute the expressions and (f79j) for N r and N into ( p3|) and 
calculate Pa . The time derivatives T and R dutifully drop out, and we get the 
relation 

-T' = R^F^APa. (80) 

When we substitute this T" back into (|7^), we can calculate P as a function of 
the canonical data: 



R'Y fP* 



Taken together, the last two equations express T" in terms of the canonical data. 
Moreover, because we know F in terms of M and R, (0), we can also determine 
the Schwarzschild mass 

M = -R^Pl - ^A~ 2 RR' 2 + -R . (82) 



Equations (|80j) - (|82j) accomplish our goal. Equation ( p2[ ) enables us to read 
the mass of the Schwarzschild black hole from the canonical data on any small 
piece of a spacelike hypersurface. It does not matter whether that piece is 
close to or far away from infinities, or even whether it lies inside or outside the 
horizon. Equations ( pp| ) and ( |3l| ) determine the difference of the Killing times 
T(ri) and Tfo) between any two points, r\ and r 2 , of the hypersurface. To 
determine T(f) itself, we need to know T at one point of the hypersurface, say, 
at the right infinity. The equations ( |80| ) - ( |82| ) are the key to our treatment of 
the Schwarzschild black holes. We shall now explore their consequences. 
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4.2 Across the Horizon 

Our assertion that equations ( ^0|) and (|8l]) determine the difference of the Killing 
times between any two points of the hypersurface requires a caveat. On the hori- 
zon, the coefficient F(r) vanishes, and the time gradient ( |8C| ) becomes infinite. 
We must show that we can propagate our knowledge of time across the horizon. 

A spacelike hypersurface must intersect both the leftgoing and the rightgoing 
branches of the horizon. Unless it passes straight through the bifurcation point, 
it has two intersections with the horizon. Inside the horizon, the hypersurface 
lies either entirely within the future dynamical region, or entirely within the 
past dynamical region. 

How does one recognize from the smooth canonical data where the hyper- 
surfaces crosses the horizon? One looks at the points where the combination 
© vanishes. The function @ can be written as a product of two factors, 

F = F + xF_, F± :=^±^. (83) 

At the horizon, at least one of the two factors must vanish. 

To find what branches of the horizon are described by the factors F±, we 
must determine the signs of several quantities at the intersection of the hyper- 
surface with the horizon. This is easily done by inspecting the Kruskal diagram. 
We chose our radial label to grow from left infinity to right infinity. Therefore, 
inside the future dynamical region T(r) increases with r, T'(r) > 0, and inside 
the past dynamical region it decreases with r, T'(r) < 0. As the hypersur- 
face is entering a dynamical region from the left static region, R(r) is falling, 
i?'(r) < 0, and as it is exiting the dynamical region into the right static region, 
R(r) is growing, R'(r) > 0. If the hypersurface passes through the bifurcation 
point, both R'(r) = and T'(r) = 0. 

In the dynamical regions, F < 0, and equation ( [sof ) tells us that T'{r) and 
Pa(^) have the same sign. Therefore, Ph(r) is positive in the future dynamical 
region, and negative in the past dynamical region. From continuity, Ph{r) > 
at intersections with the future horizon V, negative at intersections with the 
past horizon A, and zero when the hypersurface crosses the bifurcation point. 
We already know that R'(r) < when the hypersurface dives from the left 
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static region through the > part of the horizon into a dynamical region, and 
R'(r) > when it reemerges through the < part of the horizon into the right 
static region. Putting these facts together, we see that the equation F+(r) = 
defines the leftgoing branch of the horizon, and F- (r) = defines the rightgoing 
branch. When both F+(r) and F— (r) simultaneously vanish, the hypersurface 
goes through the bifurcation point. 

We can now return to the problem of determining the passage of time across 
the horizon. To be definite, let the hypersurface cross the horizon from the future 
dynamical region into the right static region at a point r : R' > 0, Pao > 0. 
Because R' > 0, we can choose R as a radial coordinate on the hypersurface 
in the vicinity of the crossing. The crossing condition F_(ro) = then implies 
(APa)q = Rq, and equation (|o|) reduces to 



f^U-i^A?^- (84) 
dR 2 R - R y ' 

for R close to Rq. From here we determine the transition of T(R) across the 
horizon: 

T(R) = -^.RoAo In \R- R \ + const. (85) 

As expected, T becomes infinite at the horizon, but in a well determined rate. 
The value of the 'constant' is determined by matching T(R) at one side of the 
horizon to its given value. 

If the hypersurface passes through the bifurcation point, T{r) changes con- 
tinuously and remains finite. We conclude that we can determine the passage 
of T(r) through the horizon from the canonical data. 

5 Mass Function and Time Gradient as Canon- 
ical Variables 

Equations ( p0[ ) - ( p^ ) were obtained from the known form of the Schwarzschild 
solution, i.e. by implicitly using the Einstein equations. Let us now forget their 
humble origin, and promote the expressions for M(r) and —T'(r) to definitions 
of two sets of dynamical variables on our phase space. Note that M(r) is a local 
functional of the canonical data, and as such it depends on r. Indeed, prior to 
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imposing the constraints and the Hamilton equations on the data, M(r) does 
not need to be constant. 

We can interpret the function M(r) as the mass content of the wormhole 
to the left of the two-sphere labeled by r. The other function, — T'(r), tells us 
the rate at which the Killing time T(r) falls with r. A remarkable feature of 
these two functions is that they form a pair of canonically conjugate variables. 
Anticipating the outcome of the proof we are going to present, we denote the 
dynamical variable —T'(r) as Pm(t). Our density notation still applies: by its 
construction (|S2"|), M(r) is a spatial scalar, while 

P M = R^F^APa (86) 

is a scalar density. 

Because the expressions for M(r) and Pm{t) do not contain Pr, they have 
vanishing Poisson brackets with R(r). Unfortunately, their Poisson brackets 
with Pr{t) do not vanish. Wc thus cannot complement the variables M(r), 
Pm{i~) by the canonical pair R{r), Pair), and get thereby a new canonical 
chart on the phase space. 

Obviously, we need to modify the momentum Pr(t) in such a way that the 
new momentum, Pr^t), will commute with M(r) and Pm (»"), but still remains 
conjugate to 

R = R . (87) 

The only way of doing this is to add to Pr(t) a dynamical variable 0(r) that 
does not depend on Pfj(r): 

Mr) =P R (r) + e(r; R,A,P A }. (88) 

To guess the correct O is tricky. Our guiding principle is that the variables 
M(r),R; Pm(t), Pr{t) should form a canonical chart whose canonical coordi- 
nates are spatial scalars, and momenta are scalar densities. This determines the 
form of the supermomcntum by the requirement that H r (r) generate DifflR. 
The same requirement had already fixed the form of the supermomentum in the 
original canonical variables. These considerations show that 

P R R' - AP A = P M M' + P R R' . (89) 
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We already know how the new canonical variables M, Pm and R depend on the 
original canonical variables. Therefore, by substituting the expressions (J8l|) - 



( p2| ) and (86) - ( p7| ) into (89) we are able to determine 6. This gives us the 
missing transformation equation for Pr : 



Pr = Pr - \r-^Pk - ^R-'F-'APa 

-R^A^F-^iAPAYiRR') - (AP A )(RR'Y) 



(90) 



By inspecting the form of the constraints, we see that Pr can be expressed as 
their linear combination: 



P R = F^tR^PAH + R'A^Hr) . 



(91) 



Our task is now clear: We must prove that the transition 

A(r),P A (r); R(r),P R (r) » M(r),P M (r); R(r),P R (r) (92) 



given by @ - @, @ - @, and © is a canonical transformation. We 
prove this by showing that the difference of the Liouville forms is an exact form: 



dr (P A (r)SA(r) + P R (r)SR(r) ) 

dr (P M (r)SM(r) + P R (r)5R(r)) = 8lo [A, Pa, R, Pr] . (93) 



The difference of the integrands is evaluated by straightforward rearrange- 
ments: 

P A SA + P R SR - P M 5M - P R SR = 



,) f AP A + ipp'ln 



RR' - AP A 



RR' + AP A 



-R8R\n 
2 



RR' + AP A 



To prove 



RR' - AP A 

|), we integrate d94) in r and argue that the boundary terms 



f 



PJPIn 



RR' + AP A 
RR' - A Pa 



(94) 



(95) 



vanish. 
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At r — ► ±00, the falloff conditions imply that A — ► 1, R — ► |r|, P' 
Pa = 0(|r|~ c ), and <5P = 0(|r|~ e ). The boundary term is of the order 

2AP A 



±1, 



P<5Pln 



1 



RR' 



5R R'- L AP A 



0(\rr) (96) 



and hence vanishes at infinities. 

There are also boundary terms at the horizon, where RSR is finite, but 
the logarithm becomes infinite. However, due to the absolute value within the 
logarithm, the infinite boundary term inside the horizon matches the infinite 
boundary term outside the horizon, and they can be considered as canceling 
each other when the integral of the derivative term is interpreted through its 
principal value. This proves (|95| ) and identifies ui : 

RR' - AP A 



oj[R, A, P A ] = / dr ( AP A + \rR' In 

J-00 V 2 



(97) 



The functional 
infinities implies 

1 



Pi?' + AP A 

?]) is well defined. The falloff of the canonical variables at 



-RR' In 
2 



RR' - AP A 



-AP a 



■0(|, 



i-(i+«h 



(98) 



RR' + AP A 

The integrand of ( |97| ) thus falls faster than | r* j 1 , which avoids the logarithmic 
singularity. Close to the horizon r = Tq, the integrand of (|97| ) behaves as 
In |r — tq I, and hence the integral from a given r to r stays finite. 

Equations ( p3| ) and (|9^) lead to the generating functional of the canonical 
transformation (p2|). The generating functional Q,[Pa,Pr; M, R] emerges when 
we introduce the old momenta and new coordinates as a new coordinate chart 
on the phase space, rewrite (p3|) in the form 



dr(A(r)SP A (r) + R(r)SR(r)) 

dr (P M (r)SM(r) + P R (r)6R(r)) = dn[P Al P R ; M, R] , (99) 



and express 



fi := 



dr (A(r)P A (r) + R(r)P R (r 
1 



w[A,P A ,P,P fi ] 



dr 



RP R + -RR' In 



RR' 


-AP A 


) 


RR' 


f AP A 





(100) 
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in the new chart. This is done by calculating A > from the mass equation 

\RR'\ 



A = 



y/R(R - 2M) + P A 2 
and substituting it back into (|100|). The result can be written in the form 



(101) 



/oo 
drR(r)P R (r 
-OO 



dr-RR' In 

2 



yR(R - 2M) + P A 2 - P A 



y/R(R - 2M) + P A 2 



Pa 



(102) 



By comparing the coefficients of the independent variations S Pa, SPr, SM, 5R in 
(|99"1) , we generate the canonical transformation (^2|) by : 



A(r) = - 6n/6P A (r) , P M (r) = - Sn/SM(r) . 
R(r) = - Sn/6P R (r) , P R (r) = - 5tt/6R(r) . 



(103) 



When resolved with respect to the new canonical variables, these equations give 
our old transformation equations (82), (|8q) , (p7|) , and (|90|). 

One should perhaps note that both integrands of ( |l02| ) fall at infinities at a 
slow rate, as 0(|r| _e ), and hence the generating functional is singular. This 
happens because we added to the well-defined functional ui an ill-defined term 



dr(P A A + P R R] 

-OO 



(104) 



The simplest way out of this difficulty is to let the momenta fall off faster than 
(3(|r|~ £ ). We have seen this can always be done for the Schwarzschild black 
hole, (||). 

From the transformation formulae Q - (||), (||) - @, @, and the weak 
falloff conditions ([l9]) - ( p2]) we easily deduce the falloff of the new canonical 
variables: 



M(t,r) 
R(t,r) 
P M {t,r) 
PR(t,r) 



= M±(t) + 0°°(|r|- 

= |r|+0°°(|rr £ ), 

= 0°°(|rr( 1+e )), 

= 00 (|r|-( 1 + £ )). 



(105) 
(106) 
(107) 
(108) 
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Recalling that Pm(t) = — T'(r), we can integrate (107) with respect to r and 
prove the earlier statement ( |55| ) about the behavior of the foliation at infinities. 

The old canonical variables are continuous (and sufficiently differentiable) 
functions of r even across the horizon. The transformation equations imply 
that the new canonical coordinates M(r) and R(r) are also continuous across 
the horizon, but this cannot be said about their conjugate momenta. Equations 
( |S6l ) and ( |9~l| ) indicate that Pm and Pr are both proportional to P" 1 . While the 
coefficients of P^ 1 are continuous, F goes to zero on the horizon, and generically 
changes its sign. As a result, Pm and Pr become infinite on the horizon, and 
generically suffer an infinite jump. 

This means that when the canonical data are such that F vanishes for some 
r (which, in particular, always happens for the Schwarzschild solution), our 
canonical transformation becomes singular. In other words, lo is not a differen- 
tiable functional of the old canonical variables. One must use the new canonical 
variables with caution. 

Except at the horizon, the canonical transformation ( p2[ ) can be inverted for 
the old variables: 

A = (F-^-FPm 2 ) 172 , (109) 
Pa - RFPm^-^-FPm 2 )- 172 , (110) 
R = R , (111) 
Pr = Pr + \p M 

+P- 1 F- 1 (F- 1 R' 2 - FP M 2 ) _1 ((RFP M )'(RR') - (RFP M )(RR')') 
+ \fPm- (112) 

In these equations, F is an abbreviation for 

F = 1-2MR -1 . (113) 

6 M(t, r) As a Constant of Motion 

Guided by the spacetime form of the Schwarzschild solution, we have introduced 
the Schwarzschild mass M(r) as a dynamical variable on our phase space. We 
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shall now prove that if the canonical data satisfy the constraints, the mass 
function M(r) does not depend on r, and if they also satisfy the Hamilton 
equations, M(r) is a constant of motion. 

Both statements follow by straightforward algebra. By differentiating the 
definition ( |S2| ) of M(r) with respect to r, we find that M'(r) is a linear combi- 
nation of the constraints (^) and (|3^): 

M' = -A~ 1 (R'H + R- 1 P A H r ) . (114) 

Because M(r) is a spatial scalar, 

{M(r),H r (r')} = M'(r) 5(r,r') , 



(115) 



equation (114) can be translated into the statement that the Poisson bracket of 
M(r) with the supermomentum weakly vanishes. 

Equations (114) and ( |9l| ) express M' and Pr as linear combinations of the 
constraints. Inversely, we can express the constraints in terms of the new canon- 
ical variables. We already know that 



H r = P R R' + P M M'. 



(116) 



The super-Hamiltonian H is then calculated from (114). By using (109) and 

F^M'R' + FP m Pr 



(110) we obtain 



H = -~, ! ^7172- ( 117 ) 

(F-!R' 2 - FP M 2 ) 

These expressions are useful for showing the closure of the Poisson bracket 



{M(r),H(r')} = -A~ 3 R' H r 5{r, r') 



(118) 



The same calculation in terms of the old variables is much more cumbersome. 
One should note that the right-hand side of (|l!8| ) does not contain the super- 
Hamiltonian, only the supermomentum. Because H{r) generates the time evo- 



lution, (11S) means that M(r) is a constant of motion. 

From equations (116) - (117) we can conclude that the Hamiltonian and 
momentum constraints 

H(r) = 0, H r (r) = (119) 
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are entirely equivalent to a new set of constraints, 



M'(r) = 0, P R (r) = 0, 



(120) 



except on the horizon. 

On the horizon, we need to be more circumspect. Let the old canonical 
variables satisfy the constraints ( |l!9| ). We know that such canonical data cor- 
respond to the Schwarzschild solution, and hence there are at most two values 
of r for which F = 0. From (114) and (|9l|) we conclude that M'{r) = every- 



where, and Pr{t) = except at the horizon points. If we insist that Pr{t) be 
a continuous function of r when the data satisfy the constraints, we conclude 
that Pr{t) — everywhere. 



Inversely, let us impose the constraints (120) on the new canonical variables. 
Except at the horizon points, equations ( |l!6| ) and ( |ll7| ) imply the old constraints 



(119). At the horizon points, F(r) = 0, and Pm(t) becomes infinite in such a 
way that F{t)Pm{t) stays finite (cf. (|8(])). We can again argue from continuity 
that the new constraints imply the old constraints even at the horizon points. 



In this sense, the constraint systems (119) and (|120|) are equivalent everywhere 



7 The Tale of Three Actions 



Written in terms of the new canonical variables M, Pm, R and Pr, the hyper- 
surface action becomes 



5 e [M,P m ,R,Pr; N,N r ] = 

dt f dr (p M (r)M(r) + P R (r)R(r) 

J — OQ 

-N(r)H(r) - N r H r {r) 



(121) 



The super-Hamiltonian H and supermomentum H r are now functions (117) 



and ( 116 ) of the new variables. The variation of (121) with respect to N and 



N r imposes the Hamiltonian and momentum constraints (HE). We have found 



that these constraints are equivalent to a new set of constraints, (120), which 



are simple functions of the new variables. The action (121) is equivalent to a 
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new action 

Sv[M,P m ,R,Pr; N m ,N r ] = 

p OO 

dt dr( P M {r)M{r) + P R (r)R(r) [ 



dt I dr (N M (r)M'{r) + N R (r)P R (r)) , (122) 

J — oo 

in which the new constraints, rather than the old ones, are adjoined to the 
Liouville form. This is done by a new set, N M (r) and N R (r), of Lagrange mul- 



tipliers. The falloff conditions (105) - (108) imply that the super-Hamiltonian 
almost coincides with M' at infinities: 

H(r) = TM'(r) + O 00 ( |r| ~^ +e > ) . (123) 

The asymptotic values of the multipliers N and N M are thus related by 

Ng (t) = TN±(t) . (124) 

The hypersurface action must again be complemented by the ADM boundary 
action (|60|). In the new variables, the ADM energy ( |63| ) is the value of the mass 
function M (r) at infinity, E± = M±. The boundary action, like the hypersurface 
action ( 122 ), is again a very simple function of the new variables: 

Sax = - J dt (N+M+ + N-M-) . (125) 

The total action is the sum 

S[M, P M , R, P R ; N M , N R ] = 5 S [M, P M , R, Pr; N m , N r ] + S ds [M; N M ] . 

(126) 

It is transparent how the boundary action cancels the boundary term obtained 
by varying M(r) in the hypersurface action. 

The lapse functions N±(t) at infinities must be treated as prescribed func- 
tions of the time parameter t. After the constraints are imposed, the boundary 
part of the total action survives as a true, ^-dependent, Hamiltonian of the black 
hole. We shall discuss this reduction process in the next section. 

The ends N±(t) of N(t, r) are freed by parametrizing the action at infinities: 



: [M- t+, r_] = - J dt (M+t+ - M_f_) . (127) 
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The total action 



S[M,P M , R,P R ; r + ,r_; N M ,N R ] = 

S s [M, P M , R, Pr; N M , N R ] + S ds [M; t+, r_] (128) 

depends now on two additional variables, T±(t), which can also be freely varied. 
It is no longer a canonical action, because r + and t_ do not come with their 
conjugate momenta. 



There are two entirely different ways in which the action ( 128) can be brought 
into canonical form. The first one is standard, the second one rather unexpected. 
Let us explain the standard method first. 



One can say that the action (128) has a mixed Hamiltonian - Lagrangian 
form, and that r±(t) is a pair of Lagrangian configuration variables. The purely 
Hamiltonian form should be reached by the Legendre dual transformation which 
complements t± by the conjugate momenta n±. However, because the action 
is linear in the velocities f±, as soon as one starts implementing this program, 
one obtains two new constraints 

C+ := n + +M + = 0, 

C- := -tt_ + M- = . 1 > 

They have vanishing Poisson brackets among themselves, and with the hyper- 
surface constraints M'(r) = = Pr(^). The resulting constraint system is thus 
first class. The new constraints C± must be adjoined to the canonical action by 
Lagrange multipliers N± : 

S[M,P m ,R,Pr] T + ,7r + ,T_,7r_; N M , N R , N+, JV_] = 

dt J dr(p M (r)M(r) + P R (r)R(r) - N M {r)M'{r) - N R {r)P R (r) ) 

dt (n + T+ + 7r_f_ - N+C+ - N-C-) . (130) 



The variation of (130) with respect to the momenta tt± leads back to equation 
(|68|). The new multipliers N± thus are what the symbols suggest: the lapse 
function at infinities. 



The price we paid for the canonical form (130) was a couple of new variables 



and a couple of new constraints. It is gratifying to learn that one can get the 
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same product for free: The mixed variables M,Pm; t+,t~ in the action (128) 
can simply be transformed into a canonical chart. It is even more gratifying 
that the new canonical variables have a clean geometric meaning: they turn out 
to be the Killing time T(r) and the mass density Pr(r) along the hypersurface, 
complemented by a canonical pair of constants of motion. 

To see how this comes about, notice that the action ( |l28[) is linear in the time 
derivatives M(r),r±. The homogeneous part of (128) thus defines a one- form 
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drP M (r)SM(r) ~ (M+6t+ - M_<5r_) 



(131) 



on (M(r),P M (r); r+,r_). 

We shall now cast ( |l3l| ) into a Liouville form. First, we replace the mass 
function M(r) by the mass at left infinity m, and by the mass density T(r): 



Inversely, 



m = M-, T(r) = M'(r). 
M(r) = m + / dr' T(r') . 



(132) 
(133) 



By introducing ( 132 ) into ( 131 ) we get 



e 



) + / drP M (r) 5m 



+ J dr ( T + 5T(r) + P M (r) / dr' 5T{r') 
+ <5(M_t_ - M+T+) . 



(134) 



To rearrange (134), we write the identity 

(f dr'P M (r')x f dr' 5T{r')\ = 

\J OO J —OO / 

dr' 5T(r') + 6T(r) / dr' P M (r') , 

- oo J OO 

which we then integrate from r ~ — ootor = oo: 



(135) 



drP M (r) / dr' 5T(r') 



dr5T(r) / dr' P M (r') . (136) 
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We immediately sec that 

6 = {i T +- T -)+ j dr'P M (r')\Sm 

+ I dr (t+- [ dr' P M {r')\ 5T(r) 

+ £(M_r_ - M+T+) . (137) 



This shows that 



to = M_, (138) 
P = (r+-r_) + / dr'P M (r') (139) 



and 



T(r) = Af'(r), (140) 
PrW = r+- / dr'P M {r') (141) 



is a canonical chart: Indeed, 

/oo 
drP r M<5r(r) (142) 
- OO 

+ 5(M_t_ - M + r+) . (143) 

differs from the Liouville form only by an exact form. 

When passing from Pm{t) = — T'(r) to T(r) we had to fix a constant of 



integration. The Killing time (141) is fixed by requiring that it matches the 
proper time r + on the parametrization clock at right infinity. By construction, 
T(r) is a spatial scalar and T(r) a scalar density. We would like to change T(r) 
into a canonical coordinate, and T(r) into a canonical momentum. This is done 
by an elementary canonical transformation (in the sense of Caratheodory |34| ) 



T(r) - P T {r)=T+- dr'P M (r'), (144) 

J OC 

P T (r) = -r(r) = -M'(r) (145) 
which sends 8 into 

/•OO 

0=p6m+ drP T (r)ST(r)+ 5lu, (146) 
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with 



/oo 
drP r (r)T(r) 
-OO 

poo pr 

M_(r_-r + )-/ dr / dr' P M (r')M'(r') . 



The last equations (146), (147) show that the transformation 
r+,r_; M(r),P M (r) i-> m,p; T(r), P T (r) 



(147) 



(148) 



given by equations (138) - (139) and (144) - ( |145| ) constructs a canonical chart 
from the originally mixed variables. 



The transformation (148) casts the parametrized action (122), (127) - (128) 
to an extremely simple canonical form 

S[m,p; T,P T ,R,P R ; N T , N R ] = 

ill ( I'm - j dr(P T {r)f(r) + P R (r)R(r) 

(149) 



dt dr(N T (r)P T (r)+N R {r)P R (r)) , 



(We gave the multiplier —N (r) a new name N (r).) In the transition from 



(128) to (149) we have discarded the total time derivative Co. Such a procedure 
does not change equations of motion, and it is used throughout classical mechan- 
ics. When applied to canonical action, it generates canonical transformations. 
Here we have used it for bringing the action to canonical form. 



Because the multipliers N T {r) and N R (r) are freely variable, the action (149) 
enforces the constraints 

P T (r) = 0, P R (r)=0. (150) 

The boundary action disappeared from ( |l49| ) and the Hamiltonian took the 
form 

/oc />oo 
drN T (r)P T (r) + / dr N R (r)P R (r) . (151) 
- oo J — oc 

It is a linear combination of constraints, and as such it weakly vanishes. 
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It is well known that the smeared super-Hamiltonian 

H[N}= drN{r)H{r) (152) 



in the Dirac-ADM action ( |3l] ) generates the change of the canonical data when 
the hypersurface is displaced by the proper time N{r) in the normal direction. 
Similarly, the smeared supermomentum 



H r [N r }= drN r (r)H r {r) (153) 

J — oo 

generates the change of the data when the hypersurface is shifted by the tan- 
gential vector N r (r). 



The Hamiltonian (151) generates a different type of displacement. The 
Hamilton equations 



T(t,r) = {T{t, r), H[N T ] } = N T {t,r), 

R(t,r) = {R(t,r), H[N T ]} = { ' 

reveal how H[N T ] displaces the hypersurface in the Kruskal diagram. It shifts 
it along the lines of constant R by the amount N T (t,r) of Killing time which 
differs from one line to another. Similarly, H[N R ] displaces the hypersurface 
along the lines of a constant T in such a way that the curvature coordinate R 



changes by the amount N R (t,r). The Hamiltonian (151) thus generates space- 
time diffeomorphisms in (T, R) . The elaboration of this general point can be 
found in Isham and Kuchaf |3f| . 

The new constraints ( |l50|) have a very simple form: a number of canonical 
momenta is set equal to zero. Locally, any system of first-class constraints 
can be brought to such a form, but it is usually impossible to find explicitly 
the necessary transformation. Our result demonstrates that this is feasible for 
Schwarzschild black holes, and that it can be done globally. The momenta which 
are required to vanish are conjugate to the embedding variables T(r) and R(r) 
which locate the hypersurface in the Kruskal diagram. 



There is a pair of canonically conjugate variables in the action (149), namely, 



m{t) and p(t), which is not subject to any constraints. However, there is no 



nonvanishing Hamiltonian in (149) which would evolve these variables. Both m 



and p are thus constants of motion. The meaning of m as the mass at infinity is 
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clear, and so is its conservation. The significance of p is at first puzzling. The 



momentum p was introduced by the transformation (13S). Due to (144), p can 
be interpreted as 

p(t)=T_(t)-r_(i). (155) 

But should not the Killing time coincide with the parametrization time also 
at left infinity, and should not thus p simply vanish? Are we not missing a 
constraint? 

The answer to this question is no . The times r± are introduced by the 
parametrization process and their origins are entirely independent. The left 
infinity does not know what the right infinity does. When we shift the origins 
by different amounts a_ = const and a + — const at the left and right infinities, 

T-{t) i-> T-(t) + a_ , r+(t) i ► T+(t) + a+ , (156) 

the parametrized action is unchanged. 

The origin of the Killing time T(r) is also arbitrary. We have chosen it 
so that T(r) matches r + at right infinity. Once defined this way, T(r) can 
be used to propagate the choice of the origin from right infinity to left infinity. 
(After the constraints are imposed, and the Schwarzschild solution is found, this 
propagation amounts to drawing the straight line across the Kruskal diagram, 
from the r + (t) = point at right infinity, through the bifurcation point of the 
horizon, and all the way up to the left infinity.) There is, however, no reason 
why the parametrization clock at left infinity should have been set to zero at 
this propagated origin. The variable p tells us the difference between the origins 
of the parametrization times at the right and left infinities. More precisely, p 
is the value of the Killing time T(r) (which is matched to the parametrization 
time at right infinity) at the origin of the parametrization time at left infinity. 
Once set, the parametrization clock and the Killing time clock run at the same 
pace, both of them measuring intervals of proper time. Therefore, it does not 
matter when we read their difference. This is the reason why p(t) of equation 
(155) is a constant of motion. 

To summarize, we have arrived at three different canonical actions describing 
the same physical system, namely, primordial black holes. The unparametrized 



canonical action (126) has a nonvanishing Hamiltonian. The two parametrized 
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canonical actions that follow have only constraints. The first of these, (130), 
has more variables and more constraints than are really necessary. Also, it does 
not disentangle the variables which are constrained to vanish from those that 



survive as true dynamical degrees of freedom. The last of our actions, (149), 
sticks to the original number of variables and constraints, and at the same time 
clearly identifies the true degrees of freedom. We believe it provides the simplest 
canonical framework for studying Schwarzschild black holes. 

8 Reduced Canonical Theory 

Each canonical action we have introduced predicts a Hamiltonian evolution. 
To compare these evolutions, we first reduce the actions to the same set of 
true degrees of freedom. The action is reduced by solving the constraints and 
substituting the solutions back into the action. 



Let us start with the unparametrized action (126). The constraint M'(r) = 



tells us that only the homogeneous mode of M (r) survives: 

M(t,r) = m(t). (157) 



By substituting (157) and Pr(V) = back into (126) we obtain 

5[m,p] = J dt (pm- (N+(t) + N-(t))mj . (158) 

The form of the reduced action enabled us to identify 

drP M {r) (159) 

as the momentum canonically conjugate to m. The reduced action has one 
degree of freedom, m, and a true time-dependent Hamiltonian 

h(t, m, p) = AT(t)m , Af(t) := N + (t) + JV_(t) . (160) 

This Hamiltonian is proportional to m, with a coefficient Af(t) which is a pre- 
scribed function of t. The Hamilton equations of motion 

rh = dh(t,m,p)/dp = 0, p = -dh(t, m, p)/<9m = -Af(t) (161) 
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indicate that m(t) is a constant of motion, but p(t) changes in time. This is 
consistent with (159) which identifies p with — (T+ — TL). The difference of the 
Killing times between the left and the right infinities stays the same only if we 
evolve the hypersurface by the lapse function which has opposite values at ±00, 
i.e. by TV = 0. 

Next, reduce the parametrized action ( 130 ) by solving both the hypersurface 
constraints and the additional constraints ( |129| ). We get 



S[m,p,T] 



dt (prh — Tm) , 



with 



T 



(162) 
(163) 



The action ( J162| ) can be obtained from the action (158) by putting 

Af(t) = T{t) . (164) 
With this replacement, the Hamilton equations of the two actions are the same, 



(161). Unlike Af(t), T(t) in action ( |162D can be varied. By varying T(t), we 
obtain once more the conservation of m(i). 

The reduction of our last action, fll49| ), by the constraints Pr{r) = = Pr{t) 
is trivial. We obtain 

S[m, P ]=J dtp rn. (165) 

We have already observed that the Hamiltonian h(m,p) vanishes, and that both 
m and p are constants of motion. 

How are the actions (158) and (162) with a nonvanishing Hamiltonian (160) 
related to the action ( 165 )? By a time-dependent canonical transformation. Let 
T(t) be a primitive function of Af(t), as in (164). Treat T(t) as a prescribed 
function of t. Take the function 



Q(t,m,p) = m(p-T(t)) 



(166) 



of the old coordinate m and the new momentum p, and let it generate a canonical 
transformation from m and p to m and p: 



p = dil(t, m,p) /dm = p — T(t) , m = d£l(t, m,p)/<9p = m . 



(167) 
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Timc-dcpcndent canonical transformations change the Hamiltonian: 



h = h + 0fi(i, m, p) jdt = T(i)m - mf(i) = . 



(168) 



Our particular generating function (166) turns our particular Hamiltonian (160) 
to zero. 



Equation (167) reproduces the definition ( |139| ) of the momentum p. The 



Hamilton equations (161) ensure that the new momentum p of equation ( 167 ) 
does not change in time. The same conclusion follows from the new Hamiltonian 



(168). The three actions generate the same dynamics. 



The canonical transformation (167) takes the canonical momentum p at t 
and transforms it to the value which it has at an instant when T(t) happens 
to vanish. The Schwarzschild mass has the same value for any t. One can thus 



view (167) as a transformation to 'initial data' 



Similarly, one can view the transition from the unparametrized action (126) 



to our final action (14E) as a time-dependent canonical transformation prior to 
the reduction. 



9 Quantum Black Holes 

The canonical action ( |149[ ) is a good starting point for the Dirac constraint quan- 
tization. The new configuration space is covered by the coordinates T(r), R(r), 
and m. The first two coordinates locate the hypersurface, the third one, m, 
is the single degree of freedom of a primordial Schwarzschild black hole. The 
constraints are 

P T (r) = 0, Pk(r)=0, (169) 

and the Hamiltonian of the system vanishes: h = 0. 

The state of the black hole on a hypersurface T(r), R(r) at the label time t 
should be described by a state functional ^(m,^ T, R] over the configuration 
space. The momenta are represented by the operators 

p= -id/dm, P T {r) = -i5/5T{r), P R {r) = -i 5 / 'SR(r) . (170) 

The Dirac rules call for imposing the constraints as operator restrictions on the 
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state functional: 



P T (r) *(m, t; T, R] = , P R (r) #(m, i; T, R] = . 



(171) 



Equations (171) imply that the state cannot depend on the embedding variables: 

* = *(m, t) . 



The state must still satisfy the Schrodinger equation 



Because h — 0, (173) ensures that W does not depend on t: 

*(m, t) = *(m) . 



(172) 



(173) 



(174) 



The state function (174) describes a superposition of primordial black holes of 
different masses. There is not much for it to do: once prepared, it stays the 
same on every hypersurface T(r), R(r) and for all t. 

Let us compare this description of states with that one which follows from 



the unparametrized canonical action (126). Let us choose the M, R] repre- 
sentation. As before, the Pr(t) — constraint implies that \& does not depend 
on R(r). The M'{r) = constraint translates into the statement that \& is an 
eigenfunction of M(r) with a constant eigenvalue M(r) = m: 

(175) 



(176) 



M'{r) * = (M(r)*J = =^ M{r) * = m * . 
In the M(r) representation 

*[M; t) = i/>(m,t)S(M(r) ~ m) , 



and we can continue working with the coefficient t/>(m, t). This coefficient must 
satisfy the Schrodinger equation with the reduced Hamiltonian (160): 

iV»(m,t) = Af(t)mip(m,t) . (177) 

Its solution is 



■0(m, t) = 4>(m) exp ( — zmT(i)) 



(178) 
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where T(t) is a primitive function to Af(t), as in (164). Unlike (174), the state 
function (178) oscillates in the T time. 

We have seen that the classical transition from the action ( |126] ) to the action 
(149) is achieved by a time-dependent canonical transformation (166) - ( |168| ). 
We want to show that the ensuing quantum theories are connected by a time- 
dependent unitary transformation 

W(i) = exp ( - iT(t)m) . (179) 

Let us start in the Heisenberg picture. The fundamental Heisenberg oper- 
ators m(£) and p(t) of the unparametrized theory depend on time, and they 
satisfy the Heisenberg equations of motion 

^ = i[m W ,h W ], ^ = i[p(*),h(C (180) 

with the Heisenberg Hamiltonian h(t) = h(t, ih(t), p(i)) . The Heisenberg states 
|0o ) refer to t = and do not depend on time. 

Change now the fundamental Heisenberg operators m(i) and p(t) into new 
fundamental Heisenberg operators rh(t) and p(t) by a general time-dependent 
unitary operator W(i) = W(i, rh(i), p(i)) : 

m(i) := W^m^W- 1 ^), (181) 
p(t) := WtfjpfiJW- 1 ^). (182) 

It is easy to show that the new Heisenberg operators form a conjugate pair, and 
that they satisfy the Heisenberg equations of motion 

d ^ = \ mm , m=i i[mm ass) 

with the new Heisenberg Hamiltonian 

- = 1 9W(t) + Tfr,^ fay ft-lus ( 184 ) 

i at 

In the Schrddinger picture, the fundamental operators rh, p and m, p become 
time- independent, while the state \<p Q ) is evolved by the respective Hamilton 
operators: 

m)) = Texp^-z^ dih(i)) |0 O ), (185) 
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i*(*): 



T exp 



-t j dth{t) ) |0 O ; 



(186) 



Here, T stands for the time ordering. We have two Schrodinger states, \ip(t) ) 
and \^f(t) ), corresponding to two alternative descriptions, m, p, h(i) and m,p,h(t), 
of the same quantum system. These states are related by 



|1r(i)) = Texp 



i J dt h(t) ) T exp 



,, I dth{t) ) 



(187) 



Apply this general scheme to our simple system. The unitary operator (179) 



yields the Heisenberg fundamental operators (|18l| ) - (182) 

m(t) = m(t) , 
p{t) = p(t)+T(t), 



(188) 
(189) 



which are related exactly as their classical counterparts (|167[) . (For simplic- 
ity, we assume that T[t = 0) = 0.) The new Heisenberg Hamilton operator 



(184) vanishes like the classical Hamiltonian (168). The Heisenberg equations 



of motion (183) then guarantee that the Heisenberg operators p(t) and m(t) are 
operator constants of motion. By ( |189| ) , the eigenvalues of the operators rh and 
rh are the same, m = m. 

The same situation can be described in the Schrodinger picture. Equation 



(187) relates the states. In m - representation, 

*(m) = *(m,t) = exp (zT(t)m) ip(m,t) 



0(m) 



(190) 



This clarifies the relation between the states (174) and (|178| 



The last of our three actions, (|130j ), has two additional constraints (129). 
The states now depend on two more configuration variables t± . The hypersur- 
face constraints reduce the states to the form ip(m, r+,r_,t). The Hamiltonian 



of the action (130) vanishes, and the Schrodinger equation implies that tp can- 
not depend on t. The reduced state function must still satisfy the (reduced) 



constraints (12E): 

C± ^(m,r + ,r_) 







( =F id/dr± + m)?/>(m, r+, r_) = . 



(191) 



These can be viewed as two Schrodinger equations in the proper times t± . Their 
solution is the state function 



ip(m, t + , T-) = 4>(m) exp ( — i m(r + — r_ 



(192) 
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This is the same state as (178), but now written in terms of the proper times r± 
rather than the label time t. Though they describe it in slightly different ways, 
our three actions lead to the same quantum dynamics. 

Because primordial black holes have only one degree of freedom m which 
is a constant of motion, their states are rather simple. Still, there are some 



interesting questions to ask. The state function ( 174 ) does not change in time. 
However, one can construct significant hypersurface-dependent operators, like 
the intrinsic and extrinsic geometry of an embedding T(r), R(r), and ask what 
their expectation values are. We defer this conceptual exercise to a later paper. 

10 Inclusion of Sources 

Schwarzschild black holes are empty vessels and what can happen to them is 
rather limited. True dynamics requires filling them with matter. This was the 
intent of the original BCMN model. Matter can propagate on the wormhole 
topology, or it can close the wormhole and change the topology of £ into 1R 3 . 
The latter case is physically more interesting. The following discussion assumes 
that the wormhole is closed. 

We must now ask whether what we have done in the vacuum can be repeated 
in the presence of matter. 

Introduce a massless scalar field propagating in the spacetime (.M,7): 

S<^; 7 ] =-{8ny 1 J d 4 X| 7 | 1/2 7 a ^,a^/3- (193) 

After the ADM decomposition and midisuperspace reduction by spherical sym- 
metry, the Lagrangian action takes the form 

S+[<f>; A,R; N,N r ] = i J dt J dr (a^ 1 AR 2 (</> - N r (/>') 2 - N A" 1 ^' 2 ) . 

° (194) 

By introducing the momentum 

7T = dL+/dj> = N~ 1 AR 2 (<j) - N r cj)') , (195) 



we cast the action ( 194 ) into canonical form by the Legendre dual transforma- 
tion: 

S*[<t>, tt; A, R; N, N r ] = f dt f dr (tt</> - NH'*' - N r H$) . (196) 



(i 
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In this process, we obtain the energy density 

^ = iA- 1 (i?-V + J RV 2 ) (197) 

and momentum density 

Hf = tt^' (198) 

of the scalar field. 

To couple the scalar field to gravity, we add the field action ( 196 ) to the 
gravitational action (|3l|) . The variation of the total action with respect to N 
and N r leads to Hamiltonian and momentum constraints on the extended phase 
space: 

H + H 4, = , H r + Hf = . (199) 

Again, up to a point transformation, this is the result obtained by BCMN pTJ. 

We arrived at the functional time formalism for the Schwarzschild black hole 
by transforming the original geometric variables into new canonical variables 
(p2[), and then into (148). Do the charms work in the presence of sources? 

The message of Section || is that (^) is a canonical transformation on the 
geometric phase space irrespective of any constraints or dynamics. Therefore, we 
can introduce the new canonical variables exactly as in the vacuum spacetime. 
(The transformation Jl48| ) needs to be modified, to accommodate the changed 
topology of E.) 

One cannot, however, expect that this immediately simplifies the constraints. 
The underlying physical reason is that the matter field curves the spacetime in 
which it propagates, and propagates in the spacetime which it curves. This is 
reflected by the presence of the metric variable A in the field energy (197) in 
the Hamiltonian constraint (199). This variable is a function (109) of the new 
canonical variables. The Hamiltonian constraint is no longer equivalent to a 
simple equation Pt{t) = —M'(r) = 0, but it provides an implicit information 
about how the energy density M'(r) depends on sources. The structure of this 
equation is presently being investigated by Dr. Joseph Romano in collaboration 
with the author. An analogous study of scalar fields coupled to a cylindrical 
gravitational wave by Braham |?6| shows that resolution of such equations is 
feasible. 
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11 Kruskal Coordinates as Phase-Space Variables 



Our main device has been the reconstruction of the curvature coordinates T and 
R from the canonical data. This turned the curvature coordinates in spacetime 
into canonical coordinates in phase space. 

Unfortunately, the Killing time T becomes infinite on the horizon, and the 
canonical transformation which leads to it has a corresponding singularity. From 
the spacetime picture one knows that the entire Schwarzschild solution can be 
covered by a single patch of spacetime coordinates, the Kruskal coordinates, 
which are well behaved on the horizon. It is natural to ask whether these coor- 
dinates, rather than the curvature coordinates, can be interpreted as canonical 
coordinates. 

The direct reconstruction of Kruskal coordinates from the canonical data is 
cumbersome. It is better to reach them via the curvature coordinates. Effec- 
tively, we are asked to reexpress the spacetime transformation (||) - ( |l3| ) as a 
point transformation on the phase space. 

Because the spacetime transformation involves exponentials, it is first nec- 
essary to turn the curvature coordinates into dimensionless quantities. This is 
done by scaling on the phase space (m, T(r), R(r); p, Pr(r), -Pr(O) • The desired 
configuration space operation mimics (g): 

T(r) - R(r) 
T(r) = -£!-, Rr = m = m. (200) 

2m 2m 

The scaled curvature coordinates T and R are dimensionless, while m keeps the 
dimension of length. 

It is important that the scaling ( |200j) be done with the Schwarzschild mass 
m at infinity rather than with the mass function M(r). Although these two 
variables coincide on the constraint surface, the coordinates T(r) and R(r) scaled 
with the mass function would not have strongly vanishing Poisson brackets, and 
hence could not be used as canonical coordinates on the phase space. 



The configuration transformation (200) can be completed into a point trans- 
formation on the phase space: 

Pt(t) = (201) 
2m 
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Inversely, 



Pk(r) = ^P, (202) 
1 r°° 

P = P-- dr(Pr(r)T(r)+P- R (r)R(r)) . (203) 



P T {r) = 2mPf(r), (204) 

P R (r) = 2mP R (r), (205) 
1 Z" 00 

p = P+— dr(P T {r)T(r)+P R (r)R(r)) . (206) 



From here, we can figure out the dimensions of the momenta. The unsealed 
coordinates T and R have the dimension of length, while the conjugate momenta 
Pt and Pr are dimensionlcss. The scaling reverts these dimensions: The scaled 
coordinates T and R are dimensionless, while the scaled momenta Pf and Pr 
have the dimension of length. Scaling does not change the dimension of the 
discrete variables: m and m, and p and p all have the dimension of length. 

Because the transformation from curvature coordinates to Kruskal coordi- 
nates is double-valued, it is better to write the transformation from Kruskal 
coordinates to curvature coordinates. The configuration part of this transfor- 
mation follows the pattern of fllO|): 

T(r) = In \V(r)\ -In \U(r)\, R(r) = n(U(r)V(r)) . (207) 



By differentiating ([L3|) with respect to T and R we obtain the Jacobi matrix 
Ut = ~\u, U R =~UF- 1 (R), (208) 
Vf = \v, V k = h/F-\R). (209) 
We took into account that 

d ^=F-^R) = {l-R-r 1 - (210) 
The completion of (|207| ) into a point transformation is straightforward: 
Pf(r) = V t f{r)Pv(r) + U f (r)Pu(r) 

= l(v(r)Py(r)-U(r)Pt,(rj), (211) 
P R (r) - V R (r)P v (r) + U R (r)Pu(r) 

= F- 1 (^(C/(r)y(r)))i(y(r)iV(r) + C/(r)P C7 (r)). (212) 
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We assume that the Kruskal variables U(r),Pu(r), V(r), Pv(r) are contin- 
uous. The transformation equations (207) and (211) - (212) reveal that the 
curvature variables R(r) and Pf( r ) will a lso be continuous, whereas T(r) and 
Pft(r) become infinite when U(r) — or V(r) = 0. 
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